Dissipative ground-state preparation of a spin chain by a structured environment 
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We propose a dissipative method to prepare the ground state of the isotropic XY spin Hamiltonian in a trans- 
verse field. Our model consists of a spin chain with nearest-neighbour interactions and an additional collective 
coupling of the spins to a damped harmonic oscillator. The latter provides an effective environment with a 
Lorentzian spectral density and can be used to drive the chain asymptotically towards its multipartite-entangled 
ground state at a rate that depends on the degree of non-Markovianity of the evolution. We also present a de- 
tailed proposal for the experimental implementation with a chain of trapped ions. The protocol does not require 
individual addressing, concatenated pulses, or multi-particle jump operators, and is capable of generating the 
desired target state in small ion chains with very high fidelities. 



I. INTRODUCTION 

The general goal of quantum information processing is to 
manipulate the information coded in a particular quantum sys- 
tem, while simultaneously trying to isolate it from undesired 
perturbations due to the external environment. For instance, 
in the fields of quantum computation [1] and quantum sim- 
ulation [2], the processing stage usually relies on a designed 
unitary evolution which, preserving the coherence of the sys- 
tem, provides a quantum state that contains the outcome of the 
computation, or the properties of the quantum phase of mat- 
ter under study. Therefore, the experimental progress in these 
fields has been typically associated to technological develop- 
ments that minimize experimental imperfections and maxi- 
mize environmental isolation [1,2]. 

However, the coupling of a system to its surrounding envi- 
ronment is not necessarily a disadvantage. In fact, the dissi- 
pation of a quantum system, if judiciously exploited, can act 
as a resource for quantum information processing. This in- 
teresting change of paradigm started with the recognition that 
dissipation can assist the generation of entanglement between 
distant atoms, either in free space [3], or trapped inside cav- 
ities [4, 5]. By measuring the presence of a photon sponta- 
neously emitted by the atoms [3, 5], or its absence [4], it is 
possible to project an initially uncorrected atomic state into a 
maximally entangled one, as has been recently experimentally 
demonstrated [6]. These schemes can be improved further by 
tailoring the atom-cavity interaction [7] or by modifying the 
measurement techniques [8], and can also be generalized to 
provide a route towards universal quantum computation [9]. 
Since the dissipation does not always render the desired re- 
sult in these approaches, the measurement is required to select 
only the successful outcomes. Therefore, the combination of 
dissipation and measurement can be considered as a proba- 
bilistic resource for quantum information. 

Another promising approach is to assess whether, by de- 
signing the system-environment coupling in a particular way, 
the dissipation could provide the desired quantum state with 
certainty. This quantum reservoir engineering finds its roots 
in the theory of laser cooling of atoms [10] and, more related 
to this work, of trapped ions [11]. Laser-cooling schemes 
try to design the system-environment coupling such that the 
dissipation yields a stationary state with reduced kinetic en- 



ergy. These ideas can be taken a step further in order to 
non-probabilistically produce stationary states that display 
non-classical aspects [12], or a certain amount of entangle- 
ment [13]. For instance, by mimicking superradiance phe- 
nomena [14], certain engineered decay channels yield partial 
entanglement in the stationary mixed state [15]. However, it 
would be highly desirable to devise dissipative protocols that 
provide maximally -entangled pure states asymptotically, the 
so-called dark states, showing fidelities comparable to those 
obtained by more standard unitary schemes [1]. 

This approach has been recently pursued by different 
groups which have shown that, provided that one can engi- 
neer a dissipation that acts quasi-locally on different parts of 
the system, it is possible to design dissipative protocols to pro- 
duce a number of paradigmatic dark states [16], or to perform 
universal dissipative quantum computation [17]. For instance, 
by engineering a particular dissipation that acts on pairs of 
adjacent atoms in an optical lattice, the system can be driven 
dissipatively into a superfluid phase [18]. In addition, as ini- 
tially proposed for Rydberg atoms [19] and first realized in ex- 
periments of trapped ions [20], by concatenating multi-qubit 
gates with a controlled dissipation of ancillary qubits, a vari- 
ety of multipartite entangled states (i.e. Bell, Greenberger- 
Home-Zeilinger, and Dicke states) have been dissipatively 
generated with reasonably high fidelities. Let us note that this 
stroboscopic time-evolution corresponds to a Markovian mas- 
ter equation in the limit of many gates and dissipation steps, 
where the dissipative jump operators correspond to multi- 
qubit jump operators. From a fundamental point of view, it 
is of interest to address whether similar dissipative protocols 
can also work in an (i) analog and (ii) global fashion (i.e. us- 
ing always-on couplings without individual addressing). Ad- 
ditionally, from a more pragmatic point of view, such global 
analog schemes would not be limited by the accumulation of 
the errors in each step of the stroboscopic protocols. However, 
finding particular schemes to provide multi-qubit jump opera- 
tors in an analog manner seems to be a tremendous task from 
both a theoretical and experimental point of view. Therefore, 
we impose a further constraint, the jump operators should be 
composed of (Hi) single-qubit operators. 

In this article, we develop an instance of a global analog dis- 
sipative protocol that generates multipartite entangled states 
corresponding to ground states of a quantum spin chain. The 



underlying idea is to engineer jump operators that are a par- 
ticular sum of single-spin operators, in analogy to the models 
of collective spontaneous emission [14]. We show that, if the 
dissipation of the spins is mediated by a common harmonic 
mode, the spin chain sees a structured Markovian environ- 
ment. By structured environment we refer to a spectral density 
exhibiting a sharp peak at a certain frequency (so called Breit- 
Wigner resonance), corresponding to a weakly damped har- 
monic oscillator. This structure will allow us to design jump 
operators in such a way that the stationary state of the system 
corresponds to the ground state of the quantum spin chain. 
In particular, we consider an isotropic XY spin chain, which 
describes a critical phase of matter. 

We describe in detail how to implement this protocol with 
trapped ions, relying on tools that have already been achieved 
experimentally. These tools are the so-called state-dependent 
forces [21-25], and sympathetic resolved- sideband cooling 
developed for quantum computation [26, 27]. We test the 
scheme for realistic parameters, and show that the fidelities 
that can be achieved are comparable to the unitary protocols 
that produce multipartite entangled states [28]. Our proce- 
dure performs well for small chains of trapped ions, which is 
actually the regime in most of the current experiments. 

This article is organized as follows. In Sec. II, we describe 
the dissipative protocol. We start by introducing the model un- 
der consideration in Sec. II A, and then move onto an analytic 
discussion, supported by numerical results, of how to engi- 
neer a structured environment that allows us to prepare mul- 
tipartite entangled states dissipatively (Sec. II B). In Sec. II C 
we explain how the method can be used to cool the system 
to the many-body ground state of the quantum spin chain. In 
Sec. Ill we show how crystals of trapped atomic ions are ide- 
ally suited for an experimental implementation of our ideas. 
Section III A describes the trapped-ion setup, a two-species 
Coulomb crystal. In Sees. IIIB, IIIC and HID we introduce 
the trapped-ion toolbox required to implement the dissipative 
protocol. Section III E contains a numerical analysis of the 
trapped-ion dissipative protocol. Finally, we present our con- 
clusions and an outlook in Sec. IV. 



II. STEADY-STATE ENTANGLEMENT OF A SPIN CHAIN 

A. Spin chain with controllable decoherence 

Let us start by introducing the model under considera- 
tion: an interacting spin-^ chain coupled to a damped bosonic 
mode (see Fig. 1(a)). This model, which shall be referred 
to as the damped spin-boson chain (DSBC) in the rest of the 
manuscript, is a many-body generalization of the system con- 
sidered recently for the dissipative generation of two-qubit en- 
tanglement [29]. The DSBC is described by the following 
master equation {h=\) 



— =^DSBC(P) 



-i[H^^H,^H,^,p]^%{p), (1) 



where p is the total density matrix, and ^dsbc is the dissipa- 
tive Liouvillian. We consider a finite chain ofN spins evolving 



under the isotropic XY Hamiltonian [30] 

N-\ 

H,= ^/(T+Oi^i+H.c. 



(2) 



where / > is an antiferromagnetic interaction strength, 
and we introduced the raising and lowering operators a^ = 
lt/)Uil = {^i~y • The bosonic Hamiltonian consists of a sin- 
gle mode 



ifb = A«<2'^a 



(3) 



where a^ and a are respectively the creation and annihilation 
operators (the oscillator frequency A^ can actually correspond 
to a detuning A^ ^ 0, as will become clear from the experi- 
mental implementation). The spins and the boson are coupled 
through 



^sb = Xl^^'^f(^ + ^^)- 



(4) 



i=l 



In the above expression, we have introduced the site- 
dependent coupling strength gi, and the Pauli matrices cr? = 
It/) (til ~ lt/)(t/l- Finally, we also consider weak damping of 
the bosonic mode via a Lindblad-type dissipator 



^b(p) = K{apa 



-a^ap) 



H.c, 



(5) 



where K* > is the damping rate. We note that all these in- 
dividual ingredients can be realized with the state-of-the-art 
technology of trapped-ion crystals, as described in Sec. III. 
The objective of this work is to understand the interplay of 
these ingredients, such that the Liouvillian ^dsbc gener- 
ates stationary entanglement between the initially uncorre- 
cted spins. So far, we can already appreciate two of the prop- 
erties of the protocol outlined in the introduction: (i) It is ana- 
log, since the couplings (2)-(5) are considered to be switched 
on during the whole protocol, (ii) It is global, since the cou- 
plings (2)-(5) address all the spins of the chain. 

Let us note that the spin-chain Hamiltonian (2) alone could 
already generate entangled states unitarily at given instants of 
time. Nonetheless, these correlations have a transient nature, 
whereas the interest of this work is the onset of stationary cor- 
relations. Therefore, we explore the interplay of this Hamil- 
tonian part with the irreversible dynamics introduced by the 
damped boson through the spin-boson coupling (4). 

In Fig. 1(a), we represent schematically the DSBC. The 
spin chain is considered to be our system y, whereas the 
bosonic mode, together with the Markovian reservoir where it 
dissipates, provide an effective structured environment (o for 
the spin chain. At first sight, the dephasing-like spin-boson 
coupling (4) seems to introduce a source of decoherence in the 
spin chain, thus hindering rather than assisting the generation 
of stationary entanglement. In the following section, we will 
show that this naive intuition is not always valid, and explain 
the subtle mechanism that allows for generation of stationary 
entanglement in a particular regime of the system parameters. 




Figure 1. Damped spin-boson chain: (a) Schematic representa- 
tion of a spin chain with an isotropic XY Hamiltonian subject to an 
additional transverse field whose intensity depends on the position 
coordinate of a damped harmonic oscillator. This oscillator mediates 
the dissipation of the spin chain, which effectively sees a structured 
environment, (b) The effect of the oscillator on the spin chain is 
described by jump operators L+ and L~ going respectively up and 
down the spectrum of spin-chain excitations with quasi-momentum 
qn and energy £n (the model parameters are defined in the text). 



B. Cooling or heating to an entangled steady state 

We consider an initial state of the system p (0) = | i/zb) ( Vb I ^^ 
Pb(0), where Pb(0) is an arbitrary state of the harmonic mode, 
and the uncorrelated spin state |v/o) = |tii2i3 '" in) has a 
single spin excitation at the left edge of the chain. In this 
section, we discuss how the system Liouvillian ^dsbc acts on 
this state, allowing for the onset of stationary entanglement. 



1. Spin-chain spectrum 

As a consequence of the XY Hamiltonian (2), the ini- 
tially localized spin excitation is exchanged between differ- 
ent neighbors. In fact, the spin chain corresponds exactly to a 
tight-binding model where the excitation hops along the sites 
of the underlying chain according to 



N-\ 



H,^i = ^sHs^s= X/|/)(/+l|+H.c.. 



(6) 



i=l 



In this expression, H^i results from the projection of the 
XY Hamiltonian onto the single-excitation subspace H^ -^ 
^s^s^s, where ^s is the projector onto an A/^-dimensional 
subspace spanned by the states | /) , such that / labels the po- 
sition of the spin excitation in the chain. This tight-binding 
model gives rise to the following band structure (see Fig. 1(b)) 

nn 



£n = 2JcOs{qn), Qn 



V--N, 



where each energy is associated to a "spin-wave excitation" 



(7) 



kn) =^^sin(^„/)|/), 



J^ = 



N^l 



(8) 



The goal of this section is to find the conditions to gen- 
erate dissipatively one of these multipartite entangled states 

TrJe^DSBc^|i/Ao)(V^o| ^Pb(O)} ^ \qn){qnl 



2. Spin-wave ladder 

We consider a modulation of the spin-boson coupling 
strength (4) in the form gi = gcos{qgi), where qg = 7l/{N-\-l). 
In the single-excitation subspace, the spin-boson coupling be- 
comes 



N-l 



/^sb,i = ^sHs^^s = Y. g{Lt^L-){a^a^), 



(9) 



n=\ 



where the ladder operators L+ = \qn){qn+\\ = {L~y are re- 
sponsible for ascending or descending along the ladder of 
spin-wave excitations (see Fig. 1(b)). Accordingly, the spin- 
boson coupling (9) connects different spin-wave excitations, 
while simultaneously creating or annihilating bosonic quanta. 
Moreover, since the harmonic mode dissipates irreversibly 
into a Markovian environment (5), the ladder operators L^ 
introduce irreversible dynamics in the spin chain. In terms of 
these ladder operators, the single-excitation spin-chain Hamil- 
tonian (6) reads 



K 



s,l 



N-l 
n=l 



lzl: 



(10) 



where for simplicity we reset the zero of energy at the bot- 
tom of the band. We show below how to generate mesoscopic 
entangled states by controlling the relative strengths of the dis- 
sipative processes that take the system up and down the spin- 
wave ladder. 



3. Irreversible dynamics in the spin chain 

Since we are interested in the stationary properties of the 
DSBC, we consider long times t ^ K~^, such that the bosonic 
mode has enough time to relax under the action of the damp- 
ing. Besides, for \g\ <C K, this relaxation is much faster than 
the process of energy exchange between the boson and the 
spins. In this limit, the boson thermalizes individually, and 
the spin chain evolves on a much slower timescale under such 
a bosonic background. 

To obtain the effect of the boson background on the slower 
spin dynamics [31], we must "integrate out" the bosonic de- 
grees of freedom from the spin-boson coupling (4). For- 
mally, this can be accomplished by projecting the density ma- 
trix of the full system into the slow spin subspace ^bP(0 = 
Pr ^Ps(0' where p^^ fulfills ^bP^ = 0, and ps{t) is the re- 
duced density matrix for the spin chain. By tracing over the 
boson degrees of freedom, one obtains 



dPs 
dt 



= ^^dTTrb{^b^sb(Oe^"^sb(^-T)^bp(0}, (11) 

in the interaction picture with respect to //s,i +ifb. where we 
have introduced ^sb(*) = — i[^sb,i ? •] and we retain up to the 
second order in g. After some algebra, and moving back to the 
Schrodinger picture, we obtain the following mediated dissi- 
pation on the spin chain 



dps 
dt 



= X(Ps) = -i[ifs,l,Ps]+^s(Ps), 



(12) 



where we have introduced the dissipation super-operator 

K 



%{ps)- 



'^coll^'^'^coll -^coir colly's 



H.c. 



(13) 



in terms of the following collective jump operators 



N-l 



N-l 



4ol = L(^t+^n), Al = '£{CLt + L-L-). (14) 

n=l n=\ 

Here, t,^ = K-/(K-±i(A„ ± A^)), being A^ = e^ - En+i < the 
energy difference between two neighboring spin- wave excita- 
tions in the spin- wave ladder (see Fig. 1(b)). 

As announced previously, the ladder operators L^ are re- 
sponsible for introducing the irreversible dynamics in the spin 
chain. In particular, they determine the collective jump opera- 
tors (14), where the adjective collective emphasizes that they 
act over all the spins in the chain. However, we remark that 
these jump operators are a sum of single-spin operators, as op- 
posed to the multi-spin nature of some other engineered dis- 
sipation protocols considered recently [16, 17, 19, 20]. With 
this discussion, we show the third property of the protocol 
outlined in the introduction: jump operators are composed 
of sums of (Hi) single-qubit operators. This draws an anal- 
ogy to the models of collective spontaneous emission [14], 
but we will show that the special dissipation mediated by the 
boson mode allows us to use them to prepare dissipatively the 

ground state of the spin chain. Let us also note that /^^ J^ con- 
tains equally-weighted dissipative processes that ascend or de- 
scend along the spin-wave ladder. Conversely, their weight 
rn = |(^^|/|(^~| can be controlled in J^J^, which will play a 
crucial role in the dissipative generation of entanglement. 




-2J l^ss) = Iqn) 



Jaiuj) A„ > 




^ss) = \qi) 




Aa<0 Jaiuj) 




-An -Ai 



Ai 



Figure 2. Tailoring the dissipative jump operators: (a,c) In the 

regime of positive oscillator detuning A,^ ~ A„, the jump operators 
that descend the ladder govern the dissipative dynamics (i.e. r„ <C 1). 
This can be understood in terms of the effective Lorentzian spectral 
density Ja of the environment seen by the spin chain, which overlaps 
minimally with the processes that climb up the ladder. (b,d) In the 
regime of negative oscillator detuning Aa ~ — A„, the jump operators 
that climb up the ladder dominate the dissipative dynamics (i.e. r„ :^ 
1), since in this case the spectral density of the environment seen by 
the spin chain overlaps minimally with the processes that descend 
the ladder. 



4. Stationary entanglement 

We are searching for the steady state pf of Eq. (12), 
^s(pI^) = 0, with three properties: being pure, displaying 
multipartite entanglement, and being unique. A pure steady 
state pf = |^ss)(^ss| necessarily belongs to the Hamiltonian 
eigenspace, which is spanned by the spin-wave excitations 
in our case |^ss) ^ span{|^„),^ = 1---N}. Additionally, the 
steady state should also belong to the kernel of the collective 

jump operator -/coiil^ss) = 0- By inspection of this jump oper- 
ator, we identify two possible regimes where this can happen, 
which depend on the boson detuning with respect to the spin- 
wave energy difference: 

K-<C+A«^A^^|A«-AJ<C|A« + AJ, r^ ^ 0, (15) 



K-<C-A«^A„^|A^ + A^|<C lA^-A^I, r^ ^ oo. (16) 



For the regime of positive detunings, we can approximate 

-^oii ^^n^n^n^ ^ jump Operator with a unique pure steady 
state corresponding to the lowest spin wave |^ss) = Wn)- 
Conversely, for negative detunings, /^^J^ ^ J^„(^+L+, and the 
steady state corresponds in this case to the highest spin-wave 



l^ss) = 1^1 )• Following [16], it can be shown that the partic- 
ular form of these limiting jump operators does not allow for 
mixed steady states. We thus get the desired result: a unique 
pure state displaying stationary multipartite entanglement. 

In Fig. 2, we represent schematically the two regimes of 
interest. For K* <C A^ ^ A„, the dissipative processes that 
climb down the spin- wave ladder dominate since | (^^ | <C | (^ ~ | . 
Effectively, the dissipation induced by the damped bosonic 
mode "cools" the spin state to the lowest-energy spin wave 
Trb{e^DSBc^|VAo)(V^o|^Pb(0)}^ |^iv)(^iv| (see Fig. 2(a)). Al- 
ternatively, for negative detunings K <C — A^ ^ A„, it is the dis- 
sipative processes that climb up the ladder which dominate, 
I ^ ~ I ^ I (^ + 1 . In this limit, the bosonic mode drives the spin 
state to the highest-energy spin wave Tr^^lQ-^^^^^^ Ixj/q) {xj/qI 
Pb(0)}^|^i)(^i| (see Fig. 2(b)). 

Let us now introduce a more physical picture that will clar- 
ify our understanding of this phenomenon. Note that the dis- 
sipation of the spin chain is mediated by the damped bosonic 
mode. Therefore, the spins do not dissipate equally at all fre- 
quencies, but more strongly at frequencies matching that of 
the bosonic mode (i.e. they do not see a totally flat environ- 
ment, but a structured one [32]). The spectral density of the 
effective environment can be recovered by realizing that the 



decay rates in Eq. (13) fulfill the following relation 
^2 






■nJaiTAn). (17) 



where the bath spectral density felt by the spin chain is 

Ja{CO) = 



;r(K-2 + (a)-A^)2)* 



(18) 



According to these expressions, the spins are under the action 
of a Lorentzian reservoir centered at the boson detuning A^ 
with a width given by the boson damping rate K. 

It is now easier to understand the two regimes in Eqs. (15)- 
(16). When conditions (15) are fulfilled, the bath spectral 
density peaks at the frequencies corresponding to the dissipa- 
tive processes descending the spin- wave ladder (see Fig. 2(c)). 
Therefore, the bath absorbs very efficiently the energy dissi- 
pated by the spin chain during its relaxation to the lowest- 
energy spin wave |^ss) = \^n)- Conversely, when the sys- 
tem parameters fulfill conditions (16), the bath spectral den- 
sity is maximal for the ascending dissipative processes (see 
Fig. 2(d)). In this case, the bath provides the required en- 
ergy to climb up the spin-wave ladder, until the chain reaches 
the highest-energy spin wave |^ss) = Wi)- We may thus con- 
clude that in both regimes, the structured reservoir singles out 
a unique spin wave as the steady state of the chain, assisting 
the generation of stationary multipartite entanglement. 

In the light of these results, we can also understand the 
regime K^ J. In this limit, the reservoir spectral density (18) 
becomes essentially flat, and both up/down processes con- 
tribute equally ^^ ^ ^~ ^ I. As a consequence, the collec- 



r(l) 



(2) 

-^coii t>ecome Hermitian, and the 



live jump operators J^^l^ 
steady state corresponds to the totally mixed state p|^ oc I. 
This regime corresponds to the naive argument of section II A, 
following which one expect that the dephasing-like term (4) 
can only decohere the spin chain. It is now clear that there 
are other regimes (15)-(16), where we can profit from a struc- 
tured environment to assist the generation of stationary entan- 
glement. 

We now comment on the time required to reach the afore- 
mentioned steady states. In the regime where the effective 
Liouvillian (12) was derived, it is given by tf ^ {\g\'^ / k)~^ , 
where \g\ <C K <^ J. For a particular experimental setup, 
where the spin couplings cannot reach arbitrarily large values, 
this preparation time may turn out to be too long for prac- 
tical purposes. We emphasize, however, that the same dis- 
sipative preparation of entangled states can be obtained in a 
non-Markovian regime where \g\ ~ K <^J, or K* <C \g\ <C /. 
Since simple master equations cannot be derived analytically, 
we shall explore this regime numerically. 




Figure 3. Approach to the asymptotic state: Time evolution of the 
error e = \ — ^, where ^ is the fidelity with the inhomogeneous 
spin-wave excitation 1^3) = (|tii) - V2|iti) + |iit))/2 for a chain 
of three spins. The curves correspond to the resonant case A = >/2/, 
while the other parameters are g = 0.01, K = 0.\ (red), g = 0.075, 
K = 0.0075 (blue), and g = K = 0.06 (green). For these numerical 
calculations, the oscillator was truncated to three levels, with trace 
preserved to order 10~^^. 



modification of the scheme also allows for the dissipative gen- 
eration of A^-partite W-like states, where the W- state is defined 
as 

|w) = ^(im---;) + im---;)+---i;;;---t)). (i9) 



Let us leave for a moment the constraint to consider purely 
global couplings, and assume that it is possible to add external 
transverse fields that act locally at the two edges of the chain 



/N-l 



H.c.) 



1 



1 



(20) 



In the subspace with one spin excitation, the ground state of 
this Hamiltonian is of the form 



\qn) 



kl*-"'"'- 



(21) 



which is locally equivalent to the above A/^-partite W-state (we 
note that the actual W-state can be obtained if the spin-spin 
coupling is ferromagnetic instead of antiferromagnetic). As 
opposed to the ground state of the spin system considered so 
far, the W-state has the property that it is fully symmetric un- 
der particle interchanges. In the state \qN), the excitation is 
equally distributed over all the sites, and each of the particles 
is equally entangled with the rest. This state can be prepared 
with the same method described before, at the price of intro- 
ducing individual addressability in the trapped-ion proposal of 
Sec. III. 



5. Dissipative preparation of W-like states 



6. Numerical analysis of the protocol 



Let us note that although the spin-wave excitations (8) are 
genuinely multipartite entangled, the weights for the excita- 
tion of each of the spins in the system are different. A small 



In this section, we analyze numerically the validity of our 
previous analytical derivations by integrating directly the mas- 
ter equation of the damped spin-boson chain (1). We obtain 



a 0.61 




C 0.6 



b 0.6 



0.5 0.2 



1 A,// 2 




d 0.6 



1 A,// 2 



Figure 4. Stationary tripartite entangled state: (a) Contour plot of the asymptotic fidelity ^oo for generating the inhomogeneous spin- 
wave excitations 1^3) = (|tii) — >/2|iti) + |iit))/2, as a function of the bosonic detuning A,^ (horizontal axis), and the spin-boson coupling 
g = K (vertical axis), for a chain of N = 3 sites (the optimal condition (15) is met for A^ = \/2/, and g <^ J), (b) The same as in (a) 
but setting a finite time tf = 10^/7. (c) Contour plot of the asymptotic fidelity ^00 for generating the homogeneous W-like states \q3) = 
(Itii) — liti) + |iit))/>/35 as a function of the bosonic detuning Aa (horizontal axis), and the spin-boson coupling g = K (vertical axis), for 
a chain ofN = 3 sites (the optimal condition (15) is met for A,^ = 7, and g <C /). (d) The same as in (c) but setting a finite time tf = 10^//. For 
these numerical calculations, the oscillator was truncated to three levels, the trace is preserved to 10~^ and the maximum fidelity is equal to 1 
within this error. In each plot, the highest fidelity is indicated by a white dot. 



the fidelity with the desired stationary entangled state 

^l^target) = I (^target |Trb{e^i^SBC^ 1 1/^0 ) (V^O | ^ Pb (0) } | ^target) | , 

(22) 

where |i//b) = |tii2i3 •••iiv), we consider that the boson 
mode is initially in the vacuum state Pb(0) = |0)(0|, and 
I ^target) IS the particular entangled state that the protocol tar- 
gets. We test the robustness of this fidelity against a non- 
optimal choice of the system parameters (15)-(16), a limited 
protocol time tf, and an increasing particle number N. 

In Figure 3, we consider the time evolution of the fidelity 
for the preparation of the lowest-energy spin- wave excitation 
of an A/^ = 3 spin chain 

\q3) = l\ni)-^\m)+l\m)- (23) 

For these calculations, we set the oscillator detuning A^ to the 
value A^ = \/2J which is resonant with the spin-wave tran- 
sitions, and vary the ratio between the spin-boson coupling 
strength g and the damping rate K. The plot indicates that, 
for similar values of the final fidelity, the convergence is much 
faster when g and K are of the same order. This implies that for 
practical purposes it is most convenient to work in the deeply 
non-Markovian regime. The optimal ratio g/ K, however, de- 
pends in general on the number of spins and the value of A^. 

It is important to note that the condition of resonance can 
be relaxed to some extent. In Figure 4 we analyze the fidelity 
as a function of the detuning A^ and of the coupling g, setting 
K = g. From our previous analysis, we expect the fidelity to be 



maximal when A^ = A^, and when g, K approach zero to ful- 
fill the requirement g^K<^J. This behavior is totally captured 
by the numerical results in Fig. 4(a), displaying the asymp- 
totic fidelity. The maximum fidelity, highlighted with a white 
dot, is indeed found in this limit and is equal to 1 within the 
numerical errors. We remark that, even when the above condi- 
tions are only partially fulfilled, the dynamics can provide the 
desired state with fidelities well above 90%. Therefore, we 
can claim that our dissipative protocol is considerably robust 
against parameter imperfections. 

Nevertheless, the time required to approach the desired tar- 
get state might be prohibitively long. Therefore, we also 
characterize the fidelity for a fixed protocol duration. In 
Fig. 4(b), we show the value of the fidelity at a finite fixed 
time tf = 10^//. As can be seen in this figure, the optimal 
choice for the oscillator detuning is still A^ = a/2/, while the 
optimal value forg=K has moved away from zero. In fact, 
the best choice now results from the interplay between the 
condition g, K* <C A^,/, and the required convergence speed. 
We emphasize that it is still possible to find parameters such 
that the achieved fidelities remain very high, ^^^ = 0.9999. 

We now address the dissipative generation of W-like states 
by adding boundary transverse fields (20). For theN = 3 spin 
chain under consideration, the target state is 

\q,) = ±{\ni)-\m)+\m))- (24) 

In Fig. 4(c)-(d), we show the asymptotic and finite-time fi- 
delities for the dissipative generation of such a multipartite 



state with homogeneously-distributed excitations. We observe 
high fidehties comparable to the inhomogeneous spin-wave 
state (8), and moreover, an analogous robustness with respect 
to a non-optimal choice of the system parameters. Let us note 
that, for A/^ = 3, there are two different transitions in the spec- 
trum of these W-type spin- waves: one is resonant at a fre- 
quency of A2 = /, and the remaining one occurs at Ai = 2/. 
Therefore, it is impossible to set an oscillator detuning at res- 
onance with both transitions simultaneously, and thus the con- 
ditions (15) are not completely fulfilled. This explains the fact 
that the W-type fidelities achieved for finite times are lower 
than those corresponding to the inhomogeneous spin waves in 
Fig. 4(b), where the two transitions have the same frequency 
Ai = A2 = V2J. 

In this section, we have analyzed numerically the validity of 
the scheme for the generation of the simplest case of multipar- 
tite entanglement, namely tripartite entangled states. A ques- 
tion that should be carefully addressed is whether the same 
scheme can be used to generate A/^-partite entangled states, and 
how large can the attained fidelities be as A^ is increased. This 
is the topic of the following section. 



7. Mesoscopic spin chains 
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We start by assessing the fulfillment of the necessary con- 
ditions (15)-(16), which rely on the energetic difference be- 
tween the dissipative processes that climb up/down the spin- 
wave ladder, as the number of spins N is increased. For 
very large spin chains N ^ 00^ the energy differences between 
neighboring spin waves |A„| < InJ /N -^ 0. This implies that 
the energetic argument selecting only processes going up or 
down the ladder can no longer hold. Indeed, limA^^cx^r^ = 
lim^^^cx. |(^ + |/|(^~| = 1, so that the protocol ceases to be op- 
erative in the thermodynamic limit of infinitely long chains. 
However, for mesoscopic spin chains, the ratio r„ can be con- 
trolled to an acceptable degree. In Fig. 5(a), we show that for 
positive detunings and N ^ 10, r^ ^ O.l for the most of the 
spin- wave excitations, whereas slightly higher values are at- 
tained for the extremal spin excitations (where |A„| oc J/N^). 
For ratios r„ on this order, we expect that the dissipative prepa- 
ration of entangled states still works with acceptable fidelities. 

To be more specific, note that the presence of different tran- 
sition frequencies A„ is generic for N >4, which represents an 
obstacle for the efficiency of the procedure. As the number of 
spins is increased, the convergence also becomes slower be- 
cause of the larger number of steps down/up the ladder to- 
wards the target state, and the lower transition frequencies 
which require a decrease in g and K. In Fig. 5(b), we show 
numerical results for the dependence of the protocol error £t^ 
achieved for a finite time tf = 10^// as a function of the num- 
ber of sites and optimizing the values of g, K and A^. We 
observe that the errors obtained for the dissipative state prepa- 
ration of the inhomogeneous spin wave Iq^) are below 10% 
for chains of up to ten sites. In comparison, the creation of 
the W-like states is worse, and such high fidelities can only be 
achieved for short chains of up to five sites. 



Figure 5. The dissipative protocol in mesoscopic chains: (a) Ratio 
fn = \^n\/\^n\ of the relative strength of the dissipative processes 
climbing up and going down the ladder as a function of quasimo- 
mentum index for different number of spins N. The values r„ > 1 
correspond to negative oscillator detunings A,^ = —^n/2 < ^ which 
select the processes up the ladder, whereas r„ < 1 corresponds to pos- 
itive oscillator detunings Aa = Aj^/2 > selecting the cooling pro- 



cesses down the ladder, (b) Error £tf = ^ 



for the dissipative 



generation of the target state in a finite time tf = 10^/7 as a function 
of the number of spins N in the chain. The blue circles represent 
the errors corresponding to the inhomogeneous spin waves |^a^) oc 
sin(^A^) \Ui ■•■;> + sin(2^A^) \iU •■•;> + ••■ + sin(A^^A^) |;U • ■ • t>. 
The red diamonds represent the error for the W-like states \qN) °^ 
\m •••;)- \m •••;) + ••• + (-1)^+^ \iU • • • t)- For the numer- 
ical calculation, the Hilbert space of the oscillator was truncated to 
four levels. 



C. Preparation of the ground state of the isotropic XY chain in 
a transverse field 



So far, our analysis has been restricted to an A/^-dimensional 
subspace of the spin chain, since the initial spin state | i/Zo) = 
Itii2i3 • • • In) contains a single excitation, and the number of 
spin excitations is preserved by the complete DSBC Liouvil- 
lian (1). In the following we explore the full 2^-dimensional 
Hilbert space of the spin chain. 



1. Jordan-Wigner jump operators 

In order to treat the full Hilbert space of the spins, we 
fermionize the spin chain via the so-called Jordan-Wigner 



transformation [33], namely 



af = 2cjc.-l, a^=cJe^^^^-^^^ = (a^)^ 



(25) 



where cj^ct are fermionic creation- annihilation operators. 
The spin-chain (2) and spin-boson (4) Hamiltonians can be 
expressed in terms of the Jordan- Wigner fermions as follows 



N-l 



ifs= I^/cJc.^i+H.c, ifsb = 1^2g,-cjc.(a + a^), (26) 



i=l 



i=l 



where we recall that the spin-boson couplings are gt = 
gcos{qgi) with qg = 7l/{N -\- 1). The next step is to express 
these Hamiltonian terms in the spin- wave basis introduced in 
Eq. (8), which leads to the following expressions 



N-l 



^s = L ^ncl^qn^ H,^ = Y, gclc^^Ja^a'^)^ll.c., (27) 
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Figure 6. Cooling to the ground state of the spin chain: (a) 

Schematic representation of the dissipative processes that cool each 
of the n^ spin excitations to the lowest non-occupied state. The sta- 
tionary state corresponds to the ground state of the isotropic XY 
Hamiltonian, with the particular filling determined by an additional 
transverse field acting as a chemical potential, (b) The same as in (a) 
but selecting the processes that climb up the ladder. The stationary 
state would correspond to the zero-temperature state of the same XY 
model in a transverse field, but with reversed coupling strengths. 



where the fermionic operators in momentum space are Cq^ = 
^L/sin(^„/)c/. In combination with the dissipative part 
given by Eq. (5), the spin-boson system is analogous to a 
damped single-mode Holstein model, a dissipative version of 
the familiar Holstein model describing electron-phonon inter- 
actions [34]. 

Concerning the single-excitation problem, we obtain the 
same formal expressions by rewriting the ladder operators in 
second-quantized form 



\qn){qn 



+ 11 



(L-y -^ LI = clc^ ^ 

\ n / n qn qn+l 



{L-J. (28) 



Accordingly, in the regime where the boson degrees of free- 
dom can be integrated out, we obtain a purely fermionic mas- 
ter equation which coincides with Eqs. (12)-(13), but with the 
collective jump operators (14) now expressed in terms of the 
Jordan-Wigner ladder operators 



-^coU 



A^-1 



= ^(L++L-), 7l^i = £(^+L+ + §-L„-). (29) 



As a result, the dissipative dynamics restricted to the single- 
excitation sector can be generalized directly to the full Hilbert 
space with arbitrary numbers of spin excitations. 



2. Effective ground-state cooling 

Let us now consider an initial state with an arbitrary number 
of spin excitations n^ <N distributed along the chain 



Due to the Pauli exclusion principle, the initial n^ excitations 
cannot all occupy the lowest energy level, with quasimomen- 
tum q^. Instead, the stationary state must be of the following 
form 



pf = |Gs)(Gs|, |Gs) 












vac). (32) 



This is precisely the ground state of the original isotropic XY 
model (2), if supplemented with a homogeneous transverse 
field ifs -^ Hs{J,jl) =Hs- (m/2) Y.i <7f • In this Hamiltonian, 
the transverse field plays the role of an effective chemical po- 
tential {1 = ^(%_(„^_i) + £N-n^)^ which is determined by the 
initial number of spin excitations. Due to the to the dissipative 
process, the excitations are distributed in the lowest available 
single-particle states (see Fig. 6(a)). 

We could also explore the steady state when the condi- 
tions (16) are fulfilled. In this limit, the jump operator 



N-l 



f^) r::,y £+^t ^ 



(33) 



n=\ 



tries to pump all the excitations to the highest-energy available 
single-particle states 



pr = |Gs)(Gs|, |Gs)=c;^---4c;|vac), (34) 

which corresponds to the ground-state of the XY Hamiltonian 
H^{—J^ —jl) with a ferromagnetic spin- spin coupling, and an 
inverted chemical potential (see Fig. 6(b)). 



\Wo) = Itit2 ••• t^si^s+ii^s+2 ■■■ in)- 



(30) 



3. Numerical results 



We would like to determine the steady-state of the spin-boson 
system if the conditions (15) are fulfilled. In this limit, the 
jump operator becomes 



. , N-l 

f^^ ^y E-c^ c 

coll L^ "^n "-qn+l^qn' 



(31) 



Let us now confirm the validity of the above argument by 
numerical analysis for small spin chains with different num- 
bers of initial excitations. In Fig. 7, we display the error in 
producing the ground-state (32) taking a fixed protocol time 
of tf = 10^//. The highest errors occur when the number 
of excitations is 1 or N —\, since the transition frequencies 




Figure 7. Scaling of the ground-state-cooling protocol: Er- 
ror £tf in producing the desired ground- state |Gs) after a fixed time 
tf = 10^/7 and for a number of spins ranging from N = 2 to N = 6. 
For each number of sites, we plot in blue circles the average over 
the different possible numbers of excitations («s ^ [^,N —\]), and 
with yellow and red squares the best and worst cases, respectively. 
For the numerical calculation, the Hilbert space of the oscillator was 
truncated to three levels. 



are lowest at these points (this is the worst case, with the er- 
rors plotted in Fig. 7 by means of red squares). With yellow 
squares, we represent the lowest attained errors. Note that the 
average errors of the protocol (blue circles, averaged over all 
the possible excitation numbers) are as low as 10~^ for chains 
up to A/^ = 6 spins, which would thus allow us to prepare dis- 
sipatively the ground state of the spin model with fidelities of 
99%. 

Let us close this section by noting that, in the general case 
with an arbitrary number of spin excitations, a different choice 
for the site-dependence of the coupling coefficients gi may 
allow for jump operators that go down/up the spin ladder in 
larger steps, and therefore favor a faster convergence to the 
asymptotic state (this appears to be specially well- suited for 
chains around the half-filling condition). However, for up to 
N = 6 spins, we have not found any sizable advantage, and the 
highest fidelities were always achieved for transitions between 
neighboring spin waves n ^ n^l. 



III. REALIZATION WITH CRYSTALS OF TRAPPED IONS 
A. Trapped-ion crystal as a spin-boson chain 

Let us describe a specific platform where the DSBC 
model (1) can be realized. We consider a system ofN-\-l ions 
confined in a linear radio-frequency trap [35], and arranged 
along a one-dimensional chain (see Fig. 8(a)). The vibra- 
tional degrees of freedom around these equilibrium positions 
can be expressed in terms of the so-called normal modes [36], 
namely 



Hn 



a,n 



'naf^na^na^ 



(35) 



where a e {x^y^z} refers to the trap axis, n e {I,. ..^N -\-l} 
labels the different modes with frequencies cOna, and bna^^lai 
are bosonic operators that annihilate/create phonons in a par- 
ticular mode. In a linear chain, the phonon branches have 
the structure shown in Fig. 8(b), such that the different 
modes spread around the trap frequencies lying in the range 
co^/2;r,^l-10MHz, cOy/27r,^l-10MHz, and coj27r ^OA- 
IMHz. A property of the longitudinal branch that will be used 
in this work is the fixed energy gap between the lowest-energy 
mode (i.e. the center-of-mass mode) and the following one. 

As shown in Fig. 8(a), N ions of the chain correspond to 
a particular atomic species with hyperfine structure, whereas 
the remaining {N -\- l)-th ion does not necessarily have hy- 
perfine splitting. For the first N ions, we select two hyper- 
fine levels of the ground-state manifold, which are referred 
to as spin states {It?)? li^)}' and have a transition frequency 
cOq in the microwave range with a negligible linewidth F ^ 
(Fig. 8(c)). For the remaining (A^+ l)-th ion, we select two in- 
ternal states {|g), 1^)} from a certain transition with frequency 
cOo, and linewidth r/2n ^10-100 kHz (Fig. 8(d)). This is the 
so-called resolved sideband limit f <C co^, which will allow for 
an efficient laser cooling close to the vibrational ground state 
for the longitudinal center-of-mass mode [26]. Let us note that 
these small linewidths can be obtained by a variety of meth- 
ods depending on the particular ion species (e.g. for weakly- 
allowed dipole transitions, for quadrupole-allowed transitions 
with an additional laser that admixes the excited state with 
that of a dipole-allowed transition, or by a Raman configura- 
tion) [37]. 

The complete degrees of freedom are thus described by the 
following master equation 

p=^{p) = -i[if,+^,+ifp,p] +#(p), (36) 

where we have introduced the Hamiltonians 



H, 



Ls^^'' 



Hn 



i=l 



2^%+l' 



(37) 



such that Tf = |t,-)(t,-| - |;,-)a,-|, and f^^^ = \e){e\ - \g){g\. 
Additionally, we have to take into account the dissipation from 
the dipole-allowed transition considering recoil effects [38]. 
For trapped ions, this can be described by the following super- 
operator #(p) = #^(p) + #^ (p), where 



n») 



(%+l •%+! 



■%+!%+! •^^.C.^ 



(38) 



where f. 



+ 



^){§\ = (%+i)^- This is the usual dissipa- 
tor of a two-level atom in an electromagnetic environment. 
The spontaneous emission also has effects on the vibrational 
degrees of freedom due to the photon recoil. To leading or- 
der in the Lamb-Dicke parameter rf^a = {G)o/c)/^/2mcOna for 
each normal mode, the jump operators create and annihilate 
phonons as given by [39]: 



-(^L+^«a)^»)^;^+l+H.C., 



(39) 
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Figure 8. Trapped-ion implementation: (a) Mixed Coulomb crystal in a linear radio-frequency trap. We consider two different types of 
ions with masses m,m, either of a different species m ^ m or a different isotope m ^ m. The crucial property is that each type has a very 
different transition frequency (Oq y^ (Oq. Colored arrows represent the combination of beams in a different configuration (i.e. R= Raman, 
tw=traveling wave, sw=standing wave) leading to the desired laser-ion interaction. Let us note that the standing-wave state-dependent forces 
can also be replaced by traveling waves without compromising severely the fidelities of the dissipative scheme, (b) Scheme of the different 
phonon branches for a two-isotope Coulomb crystal. The vibrational frequencies COna span around the trap frequency of each axis ft)«. Note 
the constant energy gap \/3ft)^ between the longitudinal center-of-mass mode, and its neighboring mode, (c) Atomic A-scheme for the N ions 
with hyperfine structure. A couple of laser beams with Rabi frequencies ^l^ 1 5 ^L^ 2 connect the two hyperfine states { 1 1) , |i) } to an auxiliary 
excited state |aux). When col^ ~ coq, the lasers drive a two-photon Raman transition, and a pair of Raman beams of this type lead to the 
state-dependent forces in Eq. (40). When cOl^ ^ (Om <C ftio, the lasers lead to a differential ac-Stark shift, such that the contribution from the 
crossed beams leads to the state-dependent force in Eq. (46). (d) Atomic scheme for the (N + l)-th ion. A driving of the transition with Rabi 
frequency ^La is red-detuned from the atomic transition cb]^^ = d)o + ^La ^ such that Al„ < 0, leading to an effective damping of the modes. 



where we have introduced f^^ °^ ^a^nai^N^inf'' ^^^ 
^N+\ n ^^^ coefficients describing the displacement of the 
A/^+ 1 ion in mode n along direction a. 

Once the vibrational and atomic degrees of freedom have 
been introduced, let us summarize the role that each of 
them will play in the realization of the damped spin-boson 
chain (1): (i) The hyperfine levels of the N ions will simulate 
the spins in the DSBC model, (ii) The 2{N ^ 1) vibrational 
modes along the x and y axes will be used to mediate the spin- 
spin interactions leading to the isotropic XY model in Eq. (2). 
(Hi) The longitudinal center-of-mass mode will play the role 
of the harmonic oscillator in the DSBC model, (iv) The atomic 
levels of the (A/^+ l)-th ion will be used to laser-cool the longi- 
tudinal center-of-mass mode. We note that the position of this 
auxiliary ion is irrelevant for the purpose of cooling, but it 
modifies the spin couplings in the chain. Moreover, it is also 
possible to use more than one auxiliary ion as an additional 
knob to control the damping rate. In the following sections, 
we will detail the different elements required for the ion-trap 
implementation of the DSBC. We emphasize that all of them 
can be realized with state-of-the-art technology. 



B. Engineering the isotropic XY spin model 

We start by addressing how to implement the isotropic XY 
model (2) in the ion chain by means of a variation of the 
M0lmer-S0rensen gate [40]. The idea is to use the so-called 
spin-dependent forces, a tool that has already been imple- 
mented in several laboratories for different purposes [21-25]. 
By combining a pair of Raman laser beams (see Fig. 8(c)), 
each of which is tuned to the red/blue first vibrational side- 
band of the hyperfine transition [41], it is possible to obtain 
the following state-dependent force 



in 



(40) 



where the "hat" indicates the interaction picture with respect 
to Hq = H(j -\- Ha ^ Hp. In the above expression, we have 
assumed that the effective wavevector of the Raman beams 
Rlc^ = kL^i — kL^2 points along the a = {x^y} axis of the 
trap (see Fig. 8(a)). Additionally, we assume that the strength 



11 



of both Raman beams is equal, such that they share a com- 
mon Rabi frequency I^l^I, and have opposite detunings, 
such that we can define a common 5na = coi^a — {<^ — ^na)- 
We have introduced the laser Lamb-Dicke parameter r{na = 
Ca • y^^i^^ / ^/2m(0yla <^ 1- The constraints on these parameters 
are I^l^ | <C (Oo^ to neglect other terms in addition to such a 
state-dependent force. Finally, we have defined the sum and 
difference of the Raman beam phases 



0a = 0L. 1 + (ka2 ' <Pa = (kaA- ^aa 



(41) 



and the PauH spin operator rf " = e-^^« 1 1/ ) {U \ + e+^^« |;,) (t,- 1 . 
In order to obtain the desired isotropic XY model, we ap- 
ply two orthogonal state-dependent forces with the following 
phases and directions 



a =x, 



0, a=y, (l)y 



K 

2' 



(42) 



Altogether, the Hamiltonian of the ion chain in the interaction 
picture becomes 

^L=I^^f,TfZ7,,e^^^e-^^-^ +^^>;TfZ7,,e^^e-^V +H.C., 

in in 

(43) 
Performing a Magnus expansion [42], the evolution turns 
out to be given by the unitary operator U{t) = exp{^eff(0}' 
where ^eff (0 displays the following form to second order 

£ieff(0 = -ifdt'HUt') - I fdt'fdt"[Hat'),Hat'% (44) 
70 2 JO Jo 

We want to derive an effective Hamiltonian from this expres- 
sion ^eff (0 ^ ~^^^s (where the superscript "ti" stands for the 
trapped-ions realization). The first-order contribution leads to 
a couple of orthogonal state-dependent displacements, which 
can only be neglected in the limit l^^^l <C 5„a [43]. In addi- 
tion, from the non-commutativity of the G^ and G^ forces, we 
get a residual spin-phonon coupling in the second-order term 
of the Magnus expansion. In order to neglect it, we have to im- 
pose a further constraint, namely \^f^{^l^y\ <C |5„x — Sjny\, 
which can be fulfilled if the trap frequencies along the x^y axes 
are sufficiently different COj^ COy. Under these constraints, an 
interacting quantum spin chain is obtained 



^s" = l4^Tj+y^TfTj, Jfj = -l^^^!^ 



jn 



)*■ 



hJ 



(45) 

Therefore, by adjusting the strengths of the Raman beams, and 
their detunings, it is possible to find a regime where Jf: = Jj-, 
and the above spin Hamiltonian corresponds to the desired 
isotropic XY model since rf Ty + t^tJ = t^t~ + ^i~^t- Let 
us remark that the emerging XY model has the peculiarity of 
displaying long-range couplings. In fact, when 0)^, (Oy ^ (0^, 
and the Raman lasers are far-detuned from the whole vibra- 
tional branch, it can be shown that these couplings decay with 
a dipolar law. The trapped-ion Hamiltonian (45) becomes a 
realization of the XY interaction (2) in the DSBC model (1) 
after the identifications 1^ ^ G^, and Jf- = Jj- ^ J. For typ- 
ical nearest-neighbor distances of zo ~l-10/im, the spin-spin 



couplings attain strengths in the J^j^^/ln ^1-10 kHz. Let 
us note that modifications of the scheme of state-dependent 
forces have been proposed to yield other quantum spin mod- 
els [43, 44], some of which have also been realized experi- 
mentally [24]. 



C. Controlling the spin-boson coupling 

Let us now address how to implement the spin-boson cou- 
pling (4) in the ion chain. The idea is again to use a A- 
beam configuration (see Fig. 8(c)), but in a different regime. 
Rather than tuning the two-photon frequencies (0\^^ to the vi- 
brational sidebands of the hyperfine transition, we impose that 
cOLc^ ~ (Onz <C (Oq. Moreover, we consider that the laser beams 
form a standing wave along the trap z-axis (see Fig. 8(a)). 
Under these constraints, the laser-ion interaction leads to a 
crossed-beam differential ac-Stark shift, which can be inter- 
preted as another state-dependent force in the G^ basis 



in 



(46) 



where the "hat" refers to the interaction picture with respect 
to Ho = Ha -\-Ha -\-Hp. Here, the parameters are defined in 
analogy to those in Eq. (40), with three important differences: 
(i) ^L^ is not the two-photon Rabi frequency of the hyperfine 
transition, but rather a differential ac-Stark shift coming from 
processes where a photon is exchanged between the pair of 
laser beams in the A configuration, (ii) The detunings are 
changed to 5nz = COl^ — ^z- (i^O Since the ion chain lies 
along the z axis, there is an additional site-dependent phase 
when shining the lasers such that Rl^ -r^ ^0. When adjust- 
ing the difference of the laser phases to be (p = 7i/2, then we 
obtain ^^ = ^r7„2^L,^^cos(kL^ -r^). Let us highlight that 
the standing- wave nature of this spin-dependent force is not 
essential for the dissipative protocol. However, it makes the 
connection with the DSBC studied in Sec. II more transparent. 
We will show in Sec. HIE that essentially the same fidelities 
can be achieved for a traveling- wave configuration, which has 
also been experimentally demonstrated [41]. 

We now exploit the particular properties of the longitudi- 
nal phonon branch (see Fig. 8(b)). In particular, we assume 
that (Ol^ is red-detuned with respect to the axial center-of- 
mass mode co^^ = coi^ — 5i^, such that 5i^ <C C0\z- Since 
the next vibrational modes are separated by a large energy 
gap \(Onz — COizl ^ V^COz^ the laser coupling to the remain- 
ing phonon branch is highly off-resonant and can be thus ne- 
glected. Accordingly, this state-dependent force gives 



^L,=l^^fiCTfZ7ue-^^i^^+H.c.. 



(47) 



We can finally move to a different picture where the above 
Hamiltonian becomes time-independent //l^ = H^ + Hl\^, 
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where 



HS 



E^^ 



U^U^U, 



^b = i:^A^f(^U + ^l)- 



(48) 



Therefore, it is clear that the center-of-mass mode plays the 
role of the harmonic oscillator {/?u,/?y ^ {a^a^} in the 
DSBC, with the laser detuning corresponding to the oscilla- 
tor detuning 5i^ ^ A^ in the DSBC model (3). In addition, 
the strength of the state-dependent force determines the spin- 
boson coupling ^:^ ^ gi in the DSBC model (4). 

Let us now comment on realistic values for the trapped- 
ion parameters. Since the laser detuning should only ful- 
fill 5ij2;r <C coijln = coj27i ?^l-10MHz, it will be easy 
to reach the required condition of 5i^/27r ~ J^j^^/ln ^l- 
10 kHz. On the other hand, we know from the previous sec- 
tion that g <C / is necessary to have an accurate pumping to 
the desired entangled state. Nevertheless, g should not be too 
small that the total preparation time becomes prohibitively 
long. A suitable choice could be g ^ 0.1-0.5 kHz. Since 
^f^ = r]u^L, cos(kL, • r^)/VN + 1, it will suffice to set the 
Rabi frequency in the ^l^ ~1-5 V^kHz. Finally, we should 
adjust the laser wavevector such that Rl^ • r^ ^ 7ri/{N-\- 1). 



D. Effective damping of the bosonic mode 

Once the implementation of the X Y model (2) and the spin- 
boson coupling (4) has been described, let us turn into the 
last required ingredient: an effective damping of the bosonic 
mode (5). As discussed previously, the idea is to exploit the 
atomic levels of the {N -\- l)-th ion to laser-cool the longitu- 
dinal center-of-mass mode. Since the resonance frequencies 
are very different, o^ 7^ (X^, the cooling lasers do not affect 
the spin dynamics of the other N ions. However, since the 
vibrational modes are collective, it is possible to sympatheti- 
cally cool the vibrations of the crystal by only acting on the 
{N^ l)-th ion. This sympathetic laser cooling [45] has been 
already realized experimentally in small crystals for quantum 
computation [26, 27]. 

Our starting point is the master equation (36) with the dis- 
sipation superoperators (38)-(39). To control the effective 
damping, we introduce a laser beam red detuned from the 
atomic transition cb^ ^ cOq -\- A^ , such that A^ < (see 
Fig. 8(d)). Let us note that the particular laser-beam con- 
figuration will depend on the particular ion species, and the 
scheme to attain the resolved-sideband limit [37]. We con- 
sider that the laser is in a traveling-wave configuration, such 
that its wavevector is aligned parallel to the trap axis Rl^ || e^ 
(see Fig. 8(a)). After expanding in series of the Lamb-Dicke 
parameter 77^^ = e^ • k^ j yjlmodnz <C 1, the laser-ion interac- 
tion leads to two different terms, the so-called carrier term 



m 



2^L, 






\(b^t 



H.C., 



(49) 



and the red and blue sideband terms 

f^l = L-^N+Unibnz + bl)a++iO-''^-^' + n.o., (50) 



~ 1 ~ ik~ -r 

where we introduced ^^+i„ = i^fjnz^t -^N+in^ ^~ '^^'' 
and the corresponding Rabi frequency Q^. ■ ^Y controlling 
appropriately these sideband terms, and their interplay with 
the atomic spontaneous emission, it is possible to tailor the 
damping of the longitudinal modes. 

Let us rearrange the full Liouvillian (36) as a sum of two 
terms ^ = ^o+-^i, where 






-i[H, + Hs+Hb + Hl,»]- 



^(•), 



-i[4-'*] 



'(•), 



(51) 



which is justified for small Lamb-Dicke parameters. To ob- 
tain the effective damping of the longitudinal center-of-mass 
mode, we use a similar formalism as in Eq. (11). However, 
in this case, the role of the fast and slow degrees of free- 
dom is interchanged. The idea is that the fastest timescale 
in the problem is given by the decay rate f of the atomic 
states of the {N -\- l)-th ion. Therefore, we must "integrate 
out" these atomic degrees of freedom, which can be accom- 
plished by projecting the density matrix of the full system into 

#p(0 = p^Vi ^P^(0. where p^^i f^lfi^^ #^p^Vi) = 0' 
and pN{t) is the reduced density matrix for the N spins and 
the chain vibrational modes. For this we use the expression 



dt 



:|^dTTr^+i{#^i(Oe^o^^i(r-T)#p(0}, 



(52) 

where the "hats" refer to the interaction picture with respect to 
Ho=H(j-\-H(j-\-Hp. After some algebra, one gets the effective 
damping of the longitudinal modes 



pN = -i[Hs -^H^^Pn] + %{Pn)^ 



(53) 



where we have introduced the bosonic dissipator 

n 

(54) 
together with the effective rates for laser cooling and heating 



K^ =Dn^Sn{±COnz)^ 



(55) 



Here, the so-called diffusion coefficient accounting for the re- 
coil heating is D^ = rnz{^N^i^N^i)ss^ and the spectral func- 
tions of the sideband terms are 

/«CX3 

Sn{(0) = j dT((/;(T)/;(0))sse^^^ F^ = #^^i^^CT++H.C., 

(56) 
which can be obtained by means of the quantum regression 
theorem. The formalism and cooling rates coincide with those 
of a single trapped ion [39], with the difference that in the 
diffusion coefficient Dn and forces Fn one must consider the 
normal vibrational modes at the position of the cooled ion. 

The idea is to work in the resolved sideband regime f <C ft)^, 
and set the laser-cooling parameters to optimize the cooling 
of the longitudinal center-of-mass mode K^ ^ K^ . In this 
regime, the cooling dominates in (54), and we are left with 
the desired damping of the bosonic mode 



m')=^iiK<z-b\Az')+<K'K-Kb 



u^L' 



-H.c. 

(57) 



13 



a 1 


B; 


■^R. 














'Bs 






0.94 








''fts 














^Q. 


»\ 












0.86 










^'e,, 

^'^u 
^^. 


0.8 










^^^E 



0.05 



0.15 0.2 




0.15 0.2 



Figure 9. Heating in the trapped-ion dissipative protocol: (a) 

Fidelity with the desired target state for a fixed time ti = 10^/7 and 
as a function of the ratio (^ = K^ /k^ between the sympathetic heat- 
ing and cooUng rates in a chain of six ions in which the central four 
represent sites of the isotropic XY spin chain. The blue dots cor- 
respond to the case of n^ = {1,3} initial spin excitations, and the 
orange dots to the case of ;is = 2 excitations. The oscillator was 
truncated to five levels, (b) Same as in (a), but considering a shorter 
time tf = 10^/7, and a traveling- wave configuration of the spin-boson 
coupling in Eq. (48). The achieved fidelities for this experimentally 
simpler configuration are similar to the ones in (a). 



E. Numerical analysis of the trapped-ion dissipative protocol 

In this last section, we will explore numerically how the 
dissipative protocol described in the part II of this manuscript 
can be implemented with the trapped-ion DSCB model in 
Eqs. (45), (48) and (57), namely 



~^ —-^DSBciP) 



-i[H^^Hf^Hi,p]^&^{p). (58) 



Therefore, the results discussed below correspond to the ac- 
tual spin- spin couplings Jf'- in an ion trap, which are not 
nearest-neighbors couplings, but display a dipolar decay with 
the cube of the distance. Furthermore, we have taken into ac- 
count that in harmonic ion traps, the interparticle distance is 
not constant over the chain. More importantly, we have con- 
sidered the effects of the always-present heating term in the 
trapped-ion effective damping (57). 

In Fig. 9(a), we show the optimal fidelity for the dissi- 
pative generation of the ideal ground state at different fill- 
ings (32), by integrating numerically the dipolar and inhomo- 
geneous trapped-ion DSBC model (58) as a function of the 
ratio between heating and cooling rate ^. We study a chain 
of A/^ + 2 = 6 ions, in which N = A have a hyperfine structure 
and play the role of the spins in the XY chain, while the two 
peripheral ions are assumed to be used for sympathetic cool- 
ing (this choice of cooling auxiliary ions at the ends of the 
chain makes the interparticle spacing slightly more homoge- 
neous). As can be seen in the figure, when ^ ^ 0, the fidelities 
with the target ground state approach 100%, which allows us 
to conclude that the differences due to the inhomogeneity of 
the chain and the dipolar range of interactions with respect to 
the ideal homogeneous and nearest-neighbor DSBC (1) do not 
compromise the fidelities severely. 

The situation is different as finite heating rates (^ > are 
considered. In Fig. 9(a), we observe that the fidelity with 
the target ground state decays essentially linear with ^. In 
the worst case considered, where ^ = 0.2 leads to a mean 
phonon number of /xi^ = 0.25, the fidelities obtained are re- 
duced to roughly ^tf ^0.8. We note that resolved sideband 
cooling of single vibrational modes has been experimentally 
achieved reaching n\z below 0.1 for single trapped ions [37]. 
Sympathetic cooling of crystals up to A^ = 4 ions has also 



Once more, after the identifications {/?u, /?{ J ^ {a, a^], and been achieved [27], reaching values as low as fii^ ^ 0.06 with 



K. ^ K, we obtain an analogous DSBC damping (5). 



Let us finally comment on the required cooling rates. Ide- 
ally, we should achieve the regime k^ g^ 0.1-0.5 kHz. Since 
such cooling is not particularly fast, it is possible to find the 
right detuning A^ < such that K^ ^0.1-0.5 kHz. Besides, 
since we are in the resolved sideband limit, the heating can be 
made much smaller, K^ <C K^ . The stationary state of a har- 
monic oscillator under the action of ^^ in Eq. (57) is a ther- 
mal state with mean number of quanta /xu = K^ /{k^ ~^\')- 
In the limit k*j^/k*j~ = ^ <C 1, such that the cooling is the dom- 
inant effect and the heating only presents a small correction, 
the equilibrium state of the center-of-mass mode is close to 
the vacuum state. 



pulsed techniques. We note that our scheme does not demand 
such an accurate ground-state cooling (see Fig. 9(a)), and we 
expect that the required mean phonon numbers ni^^OA can 
also be achieved in a continuous cooling scheme. 

Another source of experimental imperfections is given by 
slow drifts of the trap frequencies which will modify the de- 
tuning 5i^ in different experimental runs, such that the con- 
ditions (15) will not be perfectly fulfilled. Nonetheless, in 
Fig. 4 we have shown that the dissipative character of the pro- 
tocol endows it with a natural robustness with respect to non- 
optimal choices of the parameters. Therefore, one can expect 
large fidelities even in the presence of slow drifts of the trap 
frequency within the kHz-range. In any case, to alleviate im- 
perfections caused by both magnetic-field dephasing and trap 
frequency drifts, we have considered a faster protocol where 
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Figure 10. Anisotropy in the trapped-ion dissipative protocol: 

(a) Fidelity with the desired target state for a fixed time ti = 10^/7 
and as a function of the anisotropy ratio (7^ — P) / {J^ + P) for a 
chain of six ions in which the central four ions represent sites of the 
slightly-anisotropic XY spin chain. The blue dots correspond to the 
case of f7s = {I5 3} initial spin excitations, and the orange dots to the 
case of /Is = 2 excitations. The oscillator was truncated to five levels. 



the required time is ?f ~ 10^// ^1-lOms. In Fig. 9(b), we 
show that the achieved fidelities in this case are still above 
0.8 (for heating/cooling ratios below 0.2). Moreover, in this 
figure we have also considered substituting the standing-wave 
force in Eq. (48) by a traveling-wave configuration, which is 
experimentally less demanding. 

The need of the G^ state-dependent force (46) forbids the 
use of magnetic-field insensitive states (i.e. clock states). 
Therefore, the hyperfine spins will be subject to fluctua- 
tions induced by non-shielded external magnetic fields, typ- 
ically leading to dephasing in a timescale of T2 ^1-lOms. 
These timescales are much shorter than the protocol times 
tf ~ 10^// ?^l-10ms, where we have taken J /27Z ~l-10kHz. 
Nevertheless, these magnetic-field fluctuations act globally 
for standard radio-frequency traps H^ = ^AcOo{t) Y,i ^f, where 
AcOo{t) is a fluctuation of the resonance frequency due to the 
Zeeman shift (note that this might not be the case for micro- 
fabricated surface traps, where fluctuating magnetic-field gra- 
dients can also arise). Since the dynamics of the DSBC (58) 
conserves the number of spin excitations, this magnetic-field 
noise only introduces a global fluctuating phase, and thus does 
not decohere the state of the system. Therefore, as far as the 
induced spin dynamics occurs within any subspace with a con- 
served number of excitations, the dissipative protocol is robust 
to global magnetic-field noise. It is thus important to stabilize 
the intensity of the state-dependent forces to the sweet spot 



where Jfj 






and make sure that \^^A <C 5na to neglect 



any term that does not conserve the number of excitations. 
In Fig. 10, we show the achieved fidelities for anisotropics in 
Jfj/Jjj 7^ 1 . The results have been obtained numerically for the 
same regime as in Fig. 9(b), and indicate that if Jfj^Jjj differ 
by 1 part in 10^, the fidelities are not severely reduced. 



IV. CONCLUSIONS AND OUTLOOK 

We have presented a method to dissipatively generate 
multipartite-entangled states corresponding to the ground 
states of small spin chains with isotropic XY Hamiltonian in 
a transverse field. The protocol for dissipative state prepa- 
ration is interesting from a fundamental point of view, as it 
illustrates how local and even purely Markovian noise can 
assist entanglement production. Indeed, the jump operators 
required are sums of terms acting on individual sites, in con- 
trast with the few-body (quasi-local) nature of the operators 
in [16, 18]. When the noise deviates from exact Markovian- 
ity, the required time for achieving the steady state is reduced, 
which illustrates the scope of non Markovian effects for prac- 
tical purposes. We would like to emphasize that the general 
idea underlying the procedure presented in this work is not 
specific to the model considered, and might be generalized 
to other spin Hamiltonians. As an example, we introduced 
a variation that can be used to generate states locally equiv- 
alent to W- states. It would also be interesting to modify our 
dissipative protocol to prepare the ground state of a gapped 
spin model. We also believe that, for longer ion chains, the 
combination of more state-dependent forces (4) with different 
wavevectors could improve the scalability of the protocol. 

We have also explained in detail how to implement this 
method in small chains of trapped ions. In our proposal, 
two internal levels of the ions embody the spin system and 
a collective motional mode represents the damped oscilla- 
tor. The preparation procedure requires the implementation 
of spin-spin interactions using the so-called M0lmer-S0rensen 
scheme, the action of a state-dependent force, and sympathetic 
cooling, all ingredients within the capabilities of present ion- 
trap technology. Numerical simulations including different 
sources of errors indicate that the protocol can produce the 
target states with fidelities comparable to the more standard 
coherent protocols [28]. The method presented is indeed ro- 
bust against a number of experimental imperfections, and the 
implementation is simplified since only global addressing of 
the ions is required. 

Our results, however, may find a realization in different ex- 
perimental setups. For instance, the application of our ideas to 
arrays of superconducting qubits in stripline resonators seems 
feasible. Morever, since this model can also be understood as 
a quadratic fermonic model with an additional chemical po- 
tential fixing the number of particles, the results may also be 
interesting for fermionic atoms in optical lattices. 

The realization of this kind of dissipative system paves 
the way for the implementation of more complex scenarios 
to study the interplay of coherent and incoherent dynamics 
giving rise to noise-induced criticality [46]. Indeed, recent 
work in the area of non-equilibrium quantum phase transitions 
[47] has shown a number of fascinating results that could be 
demonstrated in systems of trapped ions using tools similar to 
the ones in our protocol. 
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